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Abstract. We briefly review the main physical and structural properties of Very Low-Mass 
stars. The most important improvements in the physical inputs required for the stellar mo- 
dels computations are also discussed. We show some comparisons with observational mea- 
surements concerning both the Color-Magnitude diagrams, mass-luminosity relations and 
mass-radius one, in order to disclose the level of agreement between the present theoretical 
framework and observations. 
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1. Introduction 

The interest in the study of physics at work 
in objects at the bottom of and below the 
main sequence (MS) dates back to the early 
demonstration made by Kumar (1963) that it 
has to exist a minimum total mass allowing 
a star to ignite hydrogen, the so called 'H- 
burning minimum mass' (HBMM), and that 
below this critical mass, hydrostatic equilib- 
rium is guaranteed by pressure provided by 
degenerate electrons. An additional historical 
reason for which low-mass and very low-mass 
stars (VLM) have been very important objects 
in the stellar astrophysics field is related to the 
fact that, before the results from microlens- 
ing surveys (MACHO, EROS, OGLE) were 
published, they were considered together with 
white dwarfs and sub-stellar objects, the main 
contributors to the dark matter budget in the 
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galaxy. Nowadays, there is a renewed interest 
in these stellar structures due to the fact that 
they are suitable candidates for the search of 
exoplanets. 

In order to make clearer the following dis- 
cussion, since now on we will refer with the 
term VLM stars to objects whose mass is lower 
or of the order of 0.5Mq, although in order to 
show better their peculiar structural and evolu- 
tionary properties we will make some compar- 
ison with stars in the low-mass range. 

2. The physical properties 

The thermodynamical conditions experienced 
by stellar structures in the VLM regime are 
very extreme and, possibly, among the most 
peculiar ones that stars can experience along 
their evolution. Just in order to have an idea, 
the central temperature that in the Sun is ~ 
10^ K is of the order of 4 x 10^ K for a VLM 



Cassisi: VLM stellar structures 



1 



Star with mass equal to about the HBMM value 
- i.e. ~ O.lMo at solar chemical composi- 
tion -, while the central density does increase 
from ~ lOOgcm"^ for the Sun to ~ SOOgcm"^^ 
in a O.IMq VLM object. In the same mass 
range, the temperature at the basis of the pho- 
tosphere ranges from ~ 6000 K to ~ 2500 K 
when decreasing the stellar mass, while the 
density at the same location spans a range from 
~ 10~''gcm-3 to ~ 10"^gcm-^ 

Under these thermal conditions, the cou- 
pling parameter F used for indicating how 
much the Coulomb interactions among ions 
are important is in the range between 0.1 - 
30, which means that non ideal. Coulomb in- 
teractions among ions are quite important in 
the evaluation of the Equation of State (EOS). 
When the density becomes quite large one has 
also to consider the occurrence of pressure dis- 
sociation and ionization. This condition is well 
achieved in the VLM regime and so these pro- 
cesses have to be also accounted for, as well 
as in the EOS is important to consider the pro- 
cess of formation/dissociation of molecules - 
mostly the H2 one. 

The electron degeneracy parameter t] in 
these peculiar objects is of the order of unity 
and it increases significantly when decreasing 
the stellar mass, so this implies that in the low- 
est mass tail of the VLM range there are con- 
ditions of partial electron degeneracy. 

On the basis of these evidence, we can de- 
scribe a VLM structure as an object where 
molecular H, atomic helium, and many other 
molecules (see below) are stable in the atmo- 
sphere and in the outer envelope layers, while 
the stellar interiors are mainly formed by a 
fully ionized H/He plasma; Coulomb interac- 
tions are important in the whole mass regime 
as well as the pressure ionization/dissociation 
process, while electron degeneracy increases 
with decreasing mass and becomes quite im- 
portant for a mass of about O.I2M0. It is evi- 
dent that, due to the complex thermodynamical 
conditions present in these objects, the eval- 
uation of a reliable EOS has been for long 
time a thorny problem, and only in recent 
times an accurate EOS suitable for stars in the 
VLM regime has been computed (Saumon et 
al. 1995). 



An important issue that one has to mana- 
ge when studying VLM stars is related to the 
complicated task of computing accurate model 
atmospheres which are indeed important, not 
only in order to allow a description of the 
emergent radiative flux - so aUowing to pro- 
vide model predictions about colors and mag- 
nitudes in several photometric planes -, but 
also for fixing the outer boundary conditions 
needed for solving the stellar structure equa- 
tions (Salaris & Cassisi 2005). The huge dif- 
ficulties in the VLM model atmosphere com- 
putations arise, once again from the peculiar 
physical properties which characterize these 
objects. In fact, due to the cool temperatures 
and high pressures of the outer stellar lay- 
ers, the opacity evaluations is made extremely 
complicated by the presence of a huge numbers 
of molecules whose individual contributions to 
the Rosseland mean opacity has to be properly 
evaluated. 

By simply looking to the stellar spec- 
trum of a VLM star (see Fig. 2 in AUard 
&Hauschildt (1995), and F Allard, this vol- 
ume) one easily reahze that, around and be- 
low Teff ~ 4000K, the presence of molecules 
such as H2O, CO, VO and TiO is very impor- 
tant. More in detail, TiO and VO governs the 
energy flux in the optical wavelength range, 
while H2O and CO are the dominant opacity 
source in the infrared window of the spectrum. 
Then, when the effective temperatures attain 
values lower than ~ 2500 - 2800 K, the pro- 
cess of grains condensation become extremely 
important, so increasing of orders of magni- 
tude the difficulty of a rehable opacity and, 
hence, model atmosphere computations. 

In addition, VLM stars share with white 
dwarfs an important opacitive properties: due 
to the large density and pressure in the outer 
layers the Colhsional Induced Absorption 
(CIA) on H2 molecules becomes a quite impor- 
tant contributor to the radiative flux absorption. 
Since the H2 molecule in its ground electronic 
state has no electric dipole, absorption of pho- 
tons can take place only via electric quadrupole 
transition. This is the reason for which low- 
density molecular hydrogen gas is essentially 
transparent throughout the visible and infrared 
portion of the spectrum. However, when the 
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density is large enough, each time a colh- 
sion between two particles occurs, the interact- 
ing pair such as H2 - H2, H2 - He or H2 - H, 
forms a sort of 'virtual molecule' which, be- 
cause of its nonzero electric dipole, can ab- 
sorb photons with a probability which is much 
higher than that of an isolated H2 molecule. 

When H2 CIA is efficient - as it occurs for 
stars with mass lower than about ~ O.2M0 - 
it suppresses the flux longward of 2yum and 
causes the redistribution of the emergent ra- 
diative flux towards shorter wavelengths, i.e. 
for suitable choices of the photometric bands 
VLM stars tend to appear 'bluer' with decreas- 
ing mass (Allard et al. 1997). 

As already mentioned, an additional pro- 
blem related to the computation of VLM stellar 
models is the determination of accurate outer 
boundary conditions to solve the set of inter- 
nal structure equations. It is common in stel- 
lar evolution models computations to adopt 
the boundary conditions provided by a grey 
model atmosphere - and usually this assump- 
tion provide results in quite good agreement 
with those based on more accurate and sophis- 
ticated model atmospheres (Vandenberg et al. 
2008). However, we note that the 'grey model 
atmosphere' approximation is valid only when 
all the following conditions are fulfilled: i) 
presence of an isotropic radiation field, ii) ra- 
diative equilibrium, iii) the radiative absorp- 
tion is independent on the photon frequency. 
However, in VLM stars the strong frequency- 
dependence of the molecular absorption co- 
efficients yields synthetic spectra which de- 
part severely from a frequency-averaged en- 
ergy distribution, and so the last condition is 
not clearly satisfied in these structure. 

More importantly, below ~ 5000 K, molec- 
ular hydrogen recombination in the envelope 
(H-i-H — > H2) reduces the entropy and thus the 
adiabatic gradient. This occurrence, coupled 
with the large radiative opacity and, hence, 
large value of the radiative gradient (V^d oc k), 
favors - according to the Schwarzschild cri- 
terium - the presence of a convective insta- 
bility in the atmosphere, so that convection 
penetrates deeply into the optically thin lay- 
ers (Auman (1969), Baraffe et al. (1995)), so 
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Fig. 1. The evolution of the abundance by mass of 
'He at the stellar centre, for some selected VLM 
structures and for Z = 3 x 10"'' and Y=0.23. 

radiative equilibrium is no longer satisfied^ 
AH these evidence provide compelhng argu- 
ments supporting the idea that in order to pro- 
vide a physically grounded description of the 
sub-atmospheric and envelope layers of VLM 
structures, boundary conditions provided by 
accurate, non-grey model atmospheres have 
to be adopted (Chabrier & Baraffe (2000)). 
Usually these model atmosphere predictions 
are 'attached' to the interior solution at a value 
of the optical depth r « 100, because at this 
depth the diffusive process approximation is 
fully fulfilled. However, Brocato et al. (1998) 
have shown that the effective temperature and 
the luminosity of VLM stellar models appear 
almost completely unaffected by any choice of 
T larger than unity. 

A quite important - and exclusive (when 
not accounting for the Pre-Main Sequence 
stage) - characteristic of VLM structures is that 
below a certain mass limit, they are fully con- 
vective objects. This is due to two concomitant 
processes: i) the large radiative opacity with 
characterizes both the interiors and outer lay- 

' This shortcoming is, at least, partially avoided 
by using a grey T(t) relation but performing the in- 
tegration (in t) only until an atmospheric layer sta- 
ble for convection is found. 
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Fig. 2. The evolution with time of the siuface stellar luminosity (L) in Lq, contribution of gravitation to 
the total luminosity (hp/l^f), central density (pc), central temperature (Tc), size of the convective core 
(Mcc) as fraction of the total mass, mass location of the bottom of the convective envelope (Mce), and 
central abundances by mass of H and 'He for two models at the transition between fully convective stellar 
structures: 0.4Mo (top panel) and 0.3Mq (lower panel), with Z=0.001, Y =0.23. 



ers, and that strongly increase the radiative gra- 
dient, ii) the decrease in the envelope layers of 
the adiabatic gradient due to the formation of 
molecules. The transition mass between fully 
convective stars and structure with a radiative 
core is around ~ 0.35Mo, the exact value de- 
pending on the stellar metallicity. 

It is also important to note that this con- 
vection is largely adiabatic due to the evi- 



dence that the large densities of these struc- 
tures make the heat convective transport quite 
efficient. This represents a significant advan- 
tage of VLM stars with respect more massive 
objects, because model predictions do not de- 
pend at all on the value adopted for the free 
parameter - the mixing length - entering in the 
mixing length theory for the treatment of the 
superadiabatic layers. On the other hand, the 
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fact that VLM have extended convective enve- 
lope or are fully convective has an important 
implication: the depth of the convective enve- 
lope or - for fully convective stars - the global 
thermal stratification of the structure is affected 
by the choice of the outer boundary condi- 
tions (see the discussion in Salaris & Cassisi 
(2005)). So, the choice of how to fix the bound- 
ary conditions at the bottom of the atmospheric 
layers does not aff'ect only the Teff predictions 
but also the global structural properties. 

As far as it concerns the thermonuclear 
H-burning processes important for the energy 
budget in VLM stars, due to the low central 
temperatures which characterize their interiors 
the nuclear reactions that are really important 
are: 

p-i-p^D-i-e"^-i-Ve 

p -I- D ^ ^He -I- y 

The process of destruction of ^He, 
^He(-'He,"^ He)2p is really effective only for 
T > 6 X 10*" K, this means that ■'He behaves as 
a pseudo-primary element for a large fraction 
of the whole core H-buming stage. This occur- 
rence has the important consequence that the 
definition of the 'Zero Age Main Sequence' 
for VLM stars is largely meaningless, because 
they Uve a huge fraction of their H-buming 
lifetime without attaining the equilibrium 
configuration for ^He as shown in fig. 1. The 
characteristic timescale for this process for 
stellar structures with a metallicity Z = 10"^ 
is of the order of ~ 15.8 Gyr for a 0.4Mo, 
and ~ 126 Gyr for a 0.15Mo, while for the 
same structures the core H-burning lifetime is 
~ 158 Gyr and ~ 1260 Gyr respectively, i.e. 
quite larger than the Hubble time. 

3. The structural and evolutionary 
properties 

In fig. 2 we show the time behavior of selected 
structural parameters for two VLM models 
with mass equal to 0.4 and 0.3Mo, respec- 
tively. The O.4M0 model behaves like models 
with moderately larger masses populating the 
upper portion of the MS. The increase of the 



^He abundance toward its equihbrium value in- 
creases the efficiency of the H-buming^ and the 
structure reacts decreasing both central tem- 
perature and density, which start increasing 
again only when the equilibrium value for ^He 
has been attained. However, fuUy convective 
structures as the 0.3Mo model behave quite 
differently, and central density keeps decreas- 
ing all along the major phase of H burning. 
One finds that the central density of such a 
model starts suddenly increasing again in the 
very last phases of central H burning, when the 
increased abundance of He and the correspond- 
ing decrease of radiative opacities induces a ra- 
diative shell which rapidly grows to eventually 
form a radiative core. 

Left panel of fig. 3 shows the effect of age 
on the HR diagram of a set of models. One can 
easily recognize that the age plays a negligible 
role on the H-R diagram location of stars below 
about O.5M0. It is curios that the O.3M0 model 
appears the less affected by age, less massive 
models shoving a progressively increasing sen- 
sitivity to the adopted ages. One can under- 
stand the reason for such a behaviour by con- 
sidering - as shown in fig. 1 - that the 0.3Mo 
model, at t=10Gyr, has already achieved at the 
centre the ^He equilibrium configuration, so 
that the following evolution is governed by the 
depletion of central H only (which has a very 
long timescale). On the contrary, less massive 
models keep increasing the central abundance 
of^He - an occurrence that affects the p-p chain 
efficiency - and readjusting the structure ac- 
cording to such an occurrence. A readjustment 
that in the less massive models is also ampli- 
fied by the decreasing electronic degeneracy 
induced by the decrease of central density. 

We show in the right panel of fig. 3, the 
HR diagram of VLM structures with an age of 
10 Gyr and for some selected metallicities: as 
expected as a consequence of increased opac- 
ity when the metallicity increases, the HR lo- 
cus becomes cooler and fainter with increasing 
heavy elements abundance. However, we wish 
here to emphasize the sinuous shape of the MS 
locus of VLM stars: indeed the MS loci show 



^ We remember that the p-p chain is more efficient 
when it acliieves the equilibrium. 
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Fig. 3. Left panel: The HR diagram of VLM structures for two different assumptions about the age, and 
for Z=0.0003, Y=0.23. The two arrows mark the location of the two bending points discussed in the text. 
Right panel: The HR diagram location of VLM stars for various assumptions on the metallicity. 
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Fig. 4. Left panel: The CM diagram in the optical bands for field VLM structures with known parallaxes 
(see Cassisi et al. (2000) for more details), with superimposed theoretical models for the labeled values 
of stellar metallicity. Right paneh The HST CM diagram of the MS locus of the galactic GC NGC 6397 
compared with VLM models for a suitable metallicity. 



to well-defined bending points as indicated by - the brighter point is located at 
the two arrows in the left panel of the same fig- Teff ~ 4500 K and corresponds to an 
ure: evolutionary mass of ~ O.SMq: the phys- 
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ical process at the origin of this change 
in the MS locus slope is the molecular 
hydrogen recombination. With decreasing 
total mass, at a certain point the outer stel- 
lar layers achieve the thermal conditions 
which make the H2 formation process 
quite efficient. This occurrence causes a 
decrease of the number of free particles 
in the plasma, and hence a decrease of 
the adiabatic gradient to Vad ~ 0.1 (to be 
compared with the classical value ~ 0.4 
valid for a perfect, mono-atomic gas). The 
decrease of the adiabatic gradient induced 
by H2 recombination produces a flatter 
temperature gradient and so the eff"ective 
temperature of the model is larger with 
respect a model where this process is not 
accounted for ((Copeland et al. 1970)). 
When decreasing the metallicity, this 
bending point shifts at larger Teu values, 
i.e. larger mass, because the outer layers of 
metal-poor stars are denser^, an occurrence 
that favors the formation of molecular 
hydrogen. Being the location of this point 
in the HR diagram strongly dependent 
on the thermal properties of the stellar 
envelope, its comparison with suitable 
observational constraints represents a 
formidable benchmark of the reliability of 
the adopted EOS prescriptions; 
- the faintest bending point is located at 
Teff ~ 2800 K and corresponds to an evo- 
lutionary mass of ~ O.ISM©. Its presence 
is due to the fact that decreasing the to- 
tal mass below this critical value the level 
of electron degeneracy does increase. As a 
consequence, the contribution of degener- 
ate electron pressure to the the total pres- 
sure increases, and eventually this will pro- 
duce the well-known mass-radius relation 
for fully degenerate objects, i.e. R oc M~*'^. 
On the basis of previous discussion, one 
can now predict that, when decreasing 
the metallicity, this point 'shifts' toward 
larger eff'ective temperature, i.e. larger stel- 



' Let us remember that, the lower the metallicity, 
the lower the opacity, and since dP/dr = g//<-, for a 
fixed gravity, at the same optical depth the pressure 
is larger. 



lar mass. Clearly, the location of this point 
in the HR diagram is more afl'ected by the 
opacity evaluations which control the outer 
layers pressure stratification than by the 
EOS. 

We have already defined the HBMM quan- 
tity, i.e. the minimum stellar mass that can 
attain thermal equilibrium supported by H- 
burning. Detailed numerical computations (see 
Chabrier & Baraffe (1997), Chabrier & Baraff'e 
(2000)) based on the most accurate available 
physics for VLM stars have shown that for so- 
lar metallicity the HBMM is equal to 0.075Mo. 
The value of HBMM does increase with de- 
creasing metaUicity, being equal to O.OSSMq 
for Z=0.0002. This is due - as already dis- 
cussed - to the fact that, at fixed total mass, de- 
creasing the metalhcity, the opacity decreases 
and the stellar structures become more and 
more dense so increasing the level of elec- 
tron degeneracy, a process that opposes to the 
achievement of the thermal conditions required 
for an efficient ignition of the H-buming pro- 
cesses. 

4. The observational properties 

4. 1. The Color-Magnitude diagrams 

The location of the MS loci for VLM struc- 
tures in an optical CMD is shown in fig. 4, and 
compared with some empirical measurements 
for both field dwarf stars with well known par- 
allaxes, and a galactic GC. It is worth noting 
that, while the theoretical sequences reproduce 
finely both the location and shape of the ob- 
served MS in the metal-poor regime, this is not 
true when considering solar metaUicity, field 
dwarfs. This occurrence has to be related to 
a shortcoming in the available color- r^// re- 
lation for the optical photometric bands in the 
metal-rich regime (see Cassisi et al. (2000) for 
a discussion on this issue). 

It is important to note that the morphology 
of the MS locus for VLM stars changes sig- 
nificantly when moving from the optical pho- 
tometric planes to ffie near- and far-infrared 
ones, as shown in the two panels of fig. 5. The 
most interesting evidence is that, from a few 
magnitudes below the MS Tum-ofl", the MS 
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Fig. 5. Left panel: The near-infrared CM diagram of the stellar population in a window of the galactic bulge 
(Zoccali et al. 2000) compared with model predictions for solar metallicity. Right panel: The CM diagram 
in the NICMOS HST filters FUOW and F187W for stellar models with mass equal or lower than O.SMq, 
and for selected metallicities. 



runs almost vertical in the NIR CMDs before 
suddenly shifting toward bluer colors with de- 
creasing stellar mass. This occurrence is re- 
lated to the ongoing competition - in these in- 
frared colors - between the tendency of the stel- 
lar color to become redder, due to the decreas- 
ing effective temperature and increasing radia- 
tive opacity in the optical, and the increasing 
CIA of H2 in the infrared (Saumon et al. 1994) 
which shifts back the flux to shorter wave- 
lengths. This process leads to quasi-constant 
color sequences from a mass of ~ O.SM©, cor- 
responding to Teff ~ 4500 K, for which the 
molecular hydrogen recombination starts be- 
ing eflicient to ~ 0.1 Mq. For stellar masses be- 
low this limit, the second eff'ect becomes the 
dominant one as a consequence of the larger 
densities (H2 CIA is proportional to the square 
power of the density), and this produces the 
blue loop at the very bottom of the MS. The 
blue loop becomes still more evident with de- 
creasing metallicities, due to the larger density 
which characterizes the atmosphere of metal- 
poor stellar structures. One can note that the 
same, solar composition models that fail to 
match the empirical data in the optical photo- 



metric plane, nicely reproduce the NIR-CMD 
of bulge stars. 

4.2. The mass-luminosity diagram 

Figure 6 shows the mass-luminosity relations 
for selected photometric bands ranging from 
the visual V band to the near-IR K band. 
One can easily note that the mass-luminosity 
relation becomes more and more insensitive 
to the stellar metallicity when going from 
the visual to the near-IR photometric bands, 
with the mass - Mk relation almost unaffected 
by the metallicity for mass below ~ 0.4Mo. 
This occuiTence is due to the fact that below 
Teif ~ 4000 - 4500 K, the opacity in the op- 
tical bands, dominated by TiO and VO, in- 
creases with metallicity so that the peak of 
the energy distribution is shifted toward larger 
wavelengths, in particular to the K band. This 
process causes a decreases in the flux emitted 
in the V band and conversely an increase in the 
K band with increasing metallicity. 

On the other hand, for a given mass, Tetr 
decreases with increasing metallicity, so that 
the total flux (F oc T^„) decreases. In the K 
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Fig. 6. Theoretical mass - luminosity relations in different photometric bands and for various metallicities. 
The empirical data correspond to measurements for field stars (Henry & McCarthy 1993). 



band the two effects compensate, and hence the 
mass - luminosity relation in the K band is al- 
most independent on the heavy elements abun- 
dance. However, this trends holds until the H2 
CIA does not suppress largely the flux emit- 
ted in the K band, as it occurs in the more 
metal-poor VLM stars with mass just above the 
HBMM. 

In passing, we note that the mass - magni- 
tude relations show also some inflection points, 
clearly associated to the same physical pro- 



cesses producing the bending points in the HR 
diagram. The presence of these points has to be 
properly taken into account when these mass 
- magnitude relations are used for retrieving 
the Initial Mass Function (IMF). In fact, in 
deriving the IMS the really crucial ingredient 
is the prime derivative of these relations with 
respect the stellar mass: since in correspon- 
dence of these inflection points the derivative 
shows a maximum, an accurate evaluation of 
this function is mandatory in order to avoid 
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Fig. 7. The mass - radius relation for low-mass and 
VLM structures (Cassisi et al. 2000). Empirical data 
are from Segransan et al. (2003). 

shortcoming in the IMF evaluation (Kroupa & 
Tout 1997). 

4.3. The mass - radius relation 

Thanks to the improvements in the obser- 
vational facilities, the radii of many VLM 
stars have been determined accurately. The 
study of eclipsing binaries, interferometric 
measurements with the Very Large Telescope 
Interferometer, and transit observations from 
microlensing surveys have provided a large 
sample of rehable radii and mass measure- 
ments for VLM and low-mass stars. Figure 7 
shows the comparison between some theoret- 
ical predictions and observed radii for stars 
in the range of interest: a quite good agree- 
ment seems to exist and this evidence provides 
a sound support to the reliability and accu- 
racy of last generation of stellar models for 
VLM structures. This notwithstanding, there 
are some indications (Ribas et al. 2008) that 
stellar models could really underestimate the 
stellar radii for masses below the solar one. 
Additional theoretical and observational inves- 
tigations on this topic are mandatory. 
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